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ABSTRACT 

A  number  of  important  problems  in  applied  mathematics  can  be  reduced  to  find¬ 
ing  stationary  values  of  functionals  (maxima,  minima,  and  critical  values) .  For 
functionals  defined  in  terms  of  integrals,  the  method  of  interval  integration  pro¬ 
vides  a  way  to  obtain  interval  (two-sided)  bounds  for  these  stationary  values.  As 
a  special  case  of  this  method,  upper  and  lower  bounds  for  eigenvalues  of  linear 
operators  can  be  obtained.  The  inclusion  of  stationary  values  in  intervals  is 
based  on  the  use  of  interval  functions  which  include  the  function  for  which  the 
functional  is  stationary,  and  its  derivatives.  A  simple  way  to  construct  such 
interval  functions  is  given,  and  examples  are  presented  of  a  minimum  and  an  eigen¬ 
value  problem.  The  improvement  of  initial  results  by  iteration  is  indicated. 
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SIGNIFICANCE  AND  EXPLANATION 


Solutions  to  many  problems  in  physical  mathematics  and  applied  analysis  can 
be  found  as  maxima ,  minima,  or  other  critical  points  of  certain  functionals  de¬ 
fined  in  terms  of  integrals.  Such  points  are  called  stationary  points  of  the 
functional,  and  include,  for  example,  eigenvalues  of  linear  operators.  The  solu¬ 
tion  of  these  so-called  variational  problems  by  scientific  computation  can  be  im¬ 
proved,  in  many  cases,  by  the  use  of  interval  integration,  which  yields  two-sided 
bound  (that  is,  lower  and  upper  bounds,  or  simply  interval  bounds)  for  stationary 
values.  In  order  to  apply  this  method,  one  needs  a  set  of  interval  functions 
which  enclose  a  function  for  which  the ’functional  is  stationary.  A  way  to  con¬ 
struct  such  interval  functions ,  using  the  boundary  conditions  and  reasonable 
assixnptions  on  the  highest  derivatives  is  given,  and  the  possibility  of  obtaining 
improved  bounds  by  iteration  is  indicated.  The  use  of  interval  methods  has  the 
additional  advantage  that  interval  values  can  be  assigned  to  the  boundary  condi¬ 
tions,  so  that  problems  in  which  the  boundary  conditions  sure  not  known  precisely 
can  be  studied,  or  the  response  of  a  system  to  a  range  of  conditions  can  be 
estimated. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  author  of  this  report. 


INTERVAL  BOUNDS  FOR  STATIONARY  VALUES  OF  FUNCTIONALS 

L.  B.  Rail 

1.  Variational  problems .  A  number  of  important  problems  in  physical  mathe¬ 
matics  and  applied  analysis,  particularly  the  calculus  of  variations  and  control 
theory,  reduce  to  finding  maxima  and  minima  of  functionals 

d.i)  f  *  f[y],  y  e  A, 

where  A  is  the  class  of  adml**lhle  function*  for  the  problem.  In  addition  to  the 
extremal  value * 

U.2)  I  = 

of  f,  one  may  seek  its  cAltccal  value* 

(1.3)  f*  »  f  fy*J , 

where  the  cAltlcal  point  y*  €  A  of  f  satisfies  the  Eu^CA  equation 

(1.4)  f 1 [y]  =0, 

it  being  assumed  in  this  case  that  the  G&teaux  derivative  f *  of  f  exists  on  A  [51 . 
Under  this  assumption,  extremal  point*  £,y  e  intA  such  that  £  -  f [%)  ,  7  «  f  [y] 
will  be  critical  points  of  f  [2].  For  simplicity,  extremal  and  critical  points  and 
values  of  a  functional  will  be  called  its  *tatlonaAij  points  and  values,  respective¬ 
ly.  A  vaSuatlonal  problem  for  f  on  A  is  to  find  one  or  more  of  the  pairs  (£,^)  # 
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it, y)  ,  (f*,y*),  if  such  exist. 

An  intZAvaJt  bound  for  a  stationary  value  f+  of  a  functional  f  is  simply  an 
interval  (a,bj  such  that  f*  €  (a,b],  that  is, 

(1.5)  a  5  f*  5  b. 

For  fwctionals  defined  in  terms  of  integrals,  such  as 

X1 

(1.6)  f[y]  -  /  f (x,y,y',y",...,y(n)>dx,  y  €  A, 

X0 

it  will  be  shown  that  the  method  of  interval  integration  [1] ,  [4] ,  provides  a 
way  to  obtain  interval  bounds  for  stationary  values  of  f.  In  the  variational 
problem  for  f  defined  by  (1.6) ,  the  class  A  of  admissible  functions  is  usually 
characterised  by  continuity,  differentiability,  and  boundary  conditions  on  y. 

It  should  be  noted  that  ont-iid&i  bounds  for  extremal  values  are  easy  to 
obtains  For  arbitrary  £  €  A,  one  has 

(1.7)  f  s  *l*J  5  f. 

Lower  bounds  for  minima  and  upper  bounds  for  maxima,  however,  are  often  not  easy 
to  obtain,  and  in  the  case  of  nonextremal  stationary  values  (such  as  intermediate 
eigenvalues  of  a  linear  operator) ,  one  is  often  completely  in  the  dark.  The  two- 
sided  boixids  (1.5)  furnished  by  interval  integration  are  easy  to  compute,  by  con¬ 
trast,  as  will  be  seen  below.  The  methodology  will  be  developed  for  functionals 
of  the  form  (1.6)  for  clarity,  and  its  immediate  extension  to  several  independent 
variables  will  be  presented  in  the  final  section. 

2.  Interval  integration.  Interval  analysis  [3]  is  the  branch  of  mathematics 
which  takes  real  bounded  intervals  (a,b)  as  Its  basic  units,  and  studies  trans¬ 
formations  of  them.  Its  relationship  to  real  analysis  is  somewhat  analogous  to 
that  of  complex  analysis,  since  the  reals  can  be  identified  with  the  subset  of 
intervals  which  have  equal  endpoints,  the  so-called  de^eneAO^e  intervals  x  -  [x,x] 
for  real  x.  An  intViual  function  Y  of  a  real  variable  x  assigns  the  interval 


(y.(x),y(x)I  to  each  x  in  its  interval  of  definition  X  *  The  inteAv&l 

integral  of  Y  over  X  is  the  interval 

X1  Xl_ 

(2.1)  J  Y(x)dx  -  [(LD)  /  £(x)dx  ,  <UD)  /  y(x)dx], 


where  (ID)  and  (UD)  denote  lower  and  upper  Darboux  integrals,  respectively  [1] . 
Since  these  Darboux  integrals  exist  for  all  real  functions,  it  follows  that  all 
interval  functions  are  integrable,  and  hence  integration  is  a  universal  operation 
in  interval  analysis  [1] . 

In  the  study  of  interval  transformations,  the  transformation  T  of  X  into  T(X) 
is  said  to  be  monotone  if  X  C  z  +  T(X)  C  T(Z)  ,  and  a  transformation  u  includes 
T  on  X  if  T (X)  C  u (X) ,  in  particular,  if  y  is  a  real  function,  then  the  interval 
function  Y  includes  y  on  X  if 

(2.2)  y  (X)  -  {y (x)  j  x  e  x)  C  Y (x) 

[3].  In  this  sense,  the  interval  function  Y  is  the  set  of  alt  real  functions  y 
such  that  £<x)  <  y(x)  <  y(x)  for  x  £  X,  and  one  writes  y  €  y  in  this  case.  The 
interval  integral  (2.1)  is  a  monotone  function  of  its  integrand  (1],  so  that 

xi  xi 

(2.3)  /  y (x)dx  C  /  Y (x) dx  for  y  G  Y, 


whether  or  not  y  has  a  real  (Riemann  or  Lebesgue)  integral.  The  real  interval 
of  a  real  function,  if  it  exists,  is  of  course  contained  in  its  interval  integral, 
which  always  exists  (1),  [4}.  The  calculation  of  the  interval  integral  of  Y  is 
simplified  if  the  endpoint  function*  £,y  of  Y  are  Riemann  (R)  integrable.  Then, 

(2.4)  /xY(x)dx  =  t(R)/x  ^(x)dx  ,  (R)/Xy(x)), 

so  that  inclusions  of  interval  integrals  can  be  computed  for  integrands  Z  3  Y 
with  Riemann  integrable  endpoint  functions  (11 . 
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3.  Interval  bounds.  Considering  the  integrand  of  (1.6)  to  be  a  fu..-uxon 

f (x,u„ ,u.  , . . . ,u  )  of  n  +  2  variables,  an  interval  inclusion  F (X,U  ,u  , . . . ,u  ) 

0  i  n  U  i  n 

of  it  can  be  constructed  by  interval  arithmetic  [3]  or  otherwise.  Then,  interval 
integration  provides  the  following  result. 

Theorem  3.1.  Suppose  that  £  G  A  is  a  stationary  point  of  the  functional  f 
defined  by  (1.6),  X  *  fty)  is  the  corresponding  stationary  value  of  f,  and  the 
interval  functions  YQ,  Y^ ,  . . . ,  on  X  «  txo#xi^  are  such  that  €  y^,  i  « 

0,1,..., n.  Then, 


(3.1) 


A€  [a,b]  -  /  F (x , Yg (x) ,Y^ (x) , . . . ,Y^ (x) )dx . 


This  result  provides  the  two-sided  bounds  (1.5)  for  X  immediately.  It  ap¬ 
pears,  however,  that  one  must  assisne  a  lot  about  £  and  its  derivatives  to  use 
(3.1).  In  many  cases,  one  only  has  to  assume  something  about  (for  example, 

that  it  is  bounded) ,  and  then  the  interval  functions  Y  Y, ,  Y  can  be 

n— 1  l  0 

constructed  by  use  of  interval  integration  and  the  boundary  conditions.  For  ex¬ 
ample,  suppose  that 


(3.2) 


y(n)  €  yn,  y<n-1>(*0)  e  y(n_1)  (x^  €  [c^J. 


Note  that  inteAval  boundary  conditions  can  be  prescribed.  Thus,  interval  tech¬ 
niques  can  be  useful  in  practical  problems  in  which  boundary  conditions  are  not 
known  precisely,  or  in  which  it  is  desired  to  study  the  behavior  of  a  system  over 
a  range  of  boundary  conditions. 

Indefinite  interval  integration  of  (3.2)  gives  the  functions 


x  x 

(3-3)  Yl<x)  -  [a0.B0l  +  !  Yn<t)dt,  Yr(x)  -  [Oj.B,]  ♦  /  Y  (t)dt. 

X0  X1 

Definition  3,1.  The  interval  function  Y  is  said  to  be  adnUAtZble.  for 

n 

the  boundary  conditions  (3.2)  if  [a., 8.)  €  Y  <x.  )  ,  [a„,8J  e  Y  (xj  ,  and  the 

xx  IiX  00  R0 
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intersection  Y^Cx)  nYR(x)  is  nonempty  for  x  G  X. 

Theorem  3.2.  If  the  interval  function  Y  is  admissible  for  the  boundary 

n 

conditions  (3.2)  and  9^  €  y  ,  then 

n 

(3.4)  ?<n‘1>  £  Vl  ■  W 

Proof.  By  construction,  the  interval  function  Y '  defined  by  (3.4)  con¬ 

tains  all  real  functions  g  such  that  g*  G  Y  ,  g(x  )  G  (a  ,8  ]  and  g(x  )  G  (a. , 9. J , 

no  0  0  I  II 

as  a  consequence  of  the  definition  of  the  indefinite  interval  integral  (1) .  QED. 

It  should  be  noted  that  Y  .  constructed  in  this  way  also  contains  other 

n-1 

real  functions  which  satisfy  the  boundary  conditions,  but  may  have  no  continuity 
or  differentiability  properties  at  all.  An  example  of  the  actual  construction 
of  an  interval  function  of  this  type  is  given  in  the  next  section. 

4.  The  simplest  problem  of  the  calculus  of  variations.  This  is  the  case 
n  *  1  of  (1.6)  [2] ,  and  to  simplify  matters  further,  the  boundary  conditions 

(4.1)  y(xQ)  ■=  yQ ,  yfx^  «  y1 

will  be  imposed.  The  class  A  of  admissible  functions  will  be  restricted  to  those 
for  which  y'  is  bounded,  that  is,  y*  G  [m,ra] ,  where  m,m  denote  constant  interval 
functions  with  the  corresponding  real  value.  Interval  integration  gives 


(4.2) 


yl(x)  *  yQ  +  [m,m]  (x  -  xQ)  ,  Yr(x)  *  yx  -  [m,m)  (x  -  x^  , 


and  thus  Y^ 
(4.3) 


[m,m]  is  admissible  for  (4.1)  if 

y,  - 


s  m  <  m, 


X1  ‘  X0 


The  graph  of  the  corresponding  interval  function  Y^  is  thus  a  parallelogram  with 
vertices  (xQ,y0)  and  (x^,y^),  bounded  above  by  the  intersecting  lines 

(4.4)  yl(x)  *  Y0  +  m(x  “  XQ>  »  yR(x>  a  *  5.(x  “  * 


and  below  by 


(4.5) 


x^U)  ■  yQ  *  -<x  "  xoJ '  ^r(x)  “  yi  "  m(x  “  xi) • 


Using  this  interval  function  YQ,  one  has  immediately  that 

X1 

(4.6)  f[y)  €  /  F(x,Y  (x) , [m,m] )dx 

xo 

on  the  class  A  of  functions  satisfying  (4.1)  for  which  y'  e  (m,a) .  For  example, 
suppose  that 

X1 

(4.7)  fly)  -  /  /I  +  (y’)*dx, 

X0 

and  one  seeks  X  *  min  f[y]  on  A.  Since 

(4.0)  (y')2  €  [o,max{m2 ,m2}) , 

which  gives 

xi _  _ _ 

(4.9)  X  €  [Xx  -  XQ  ,  d],  d  -  /  /l  +  mTdx  -  /(x^  -  xQ)  z  +  (y^  -  yQ)'* , 

X0 

since  each  admissible  contains  y*  *  m,  and  thus  each  VQ  contains  the  degen¬ 
erate  interval  (real)  function 

(4.10)  yB(x)  -  yQ  ♦  m(x  -  xQ)  , 

for  which  the  value  Fty^]  *  d  is  attained,  so  that  X  s  d.  The  lower  bound  in 
(4.9) ,  X  »  x^  -  xQ,  can  be  attained  by  a  step  function  s  €  YQ  if  0  €  Y^.  in 
variational  problems,  additional  restrictions  are  often  imposed  on  the  elements 
of  Yq  (continuity,  boundedness  of  derivatives)  to  eliminate  solutions  such  as 
this,  and  one  works  with  the  fractions  in  YQ  n  A  121 j  in  this  class,  however, 
two-sided  bounds  may  not  be  available  as  they  are  in  interval  analysis. 

5.  Eigenvalue  problems.  For  self ad joint  linear  operators  A  in  a  Hilbert 
space  H,  its  eigenvalue*  are  critical  values  of  the  Rayleigh  quoliznt 

(5.1)  R(y)  -  (Ay,y)/(y,y) ,  y  ?  0. 

An  eigenvalue  X  satisfies  the  Euler  equation 

-  6  - 


(5.2) 


Ay  -  Xy  *  o,  y  ¥  0 


and  the  corresponding  critical  points  y  in  a  function  space  H  are  called 
function*  of  A  belonging  to  X  (2).  If  the  inner  product  (  ,  )  in  H  is  defined 
in  terms  of  integrals,  then  interval  integration  can  be  applied  as  above  to  find 
lower  and  upper  bounds  for  eigenvalues  of  A.  If  Y^,  Y^  are  interval  functions 
such  that  ?  €  y  ,  0  €  Y  #  and  A?  €  for  seme  eigenfunction  £  of  A,  then 

(5.3)  X  e  (YlfY0)/(Y0#Y0)  =  [a,bl 

for  the  corresponding  eigenvalue  X,  thus  giving  an  interval  bound.  Once  again, 
the  interval  functions  YQ,  Y^  are  to  be  determined  in  some  way,  perhaps  on  the 
basis  of  an  approximate  solution  of  (5.2).  If  A  is  an  integral  operator,  then 
one  may  use  interval  integration  to  get  Y^  *  AYq.  0,1  t*le  ot^er  hand,  if  A  is  a 
differential  operator,  it  may  be  possible  to  obtain  Yq  from  Y^  by  use  of  interval 
integration  and  the  boundary  conditions,  as  before.  For  example,  suppose  that 

(5.4)  Ay  =  -y",  y (0)  »  y ( it)  *  o. 

Since  eigenfunctions  are  determined  by  (5.1)  and  (5.2)  only  up  to  a  multiplicative 
constant,  it  is  useful  to  introduce  a  normalization  condition  which  excludes  y  *  0 


(5.8) 


*0 


(x)  * 


x<-  -  — )  ,  0  <  x  <  */2, 

4 


(TT  -  x)  (“  -  —  "  )  ,  7T/2  £  x  <  TT. 

4  3tt 


Computation  with  these  interval  functions  gives 


(5.9) 


(vv  -t'H1'  (yo’V  =  11 ' 


and  thus  eigenvalues  A  belonging  to  eigenfunctions  y  G  Y^  will  be  elements  of  the 
intwoat  Rayleigh  quotient  <y1#y0)/(Yo#y  )#  that  is, 

(5.10)  X  €  A0  "  2^5  *  I7]  C  l0*4052  '  2.5033]  . 


In  this  case,  y(x)  =  sin  x  is  the  only  eigenfunction  of  A  contained  in  YQ,  and 

(5.10)  provides  lower  and  upper  bounds  for  the  corresponding  eigenvalue  A  *  1. 

The  endpoint  functions  (5.5)  are  crude  approximations  to  sinx,  and  it  can  be 

noted  that  (5.7)  and  (5.8)  define  an  interval  function  which,  when  normalized, 

is  smaller  than  and  bounded  by  better  approximations  to  the  eigenfunction.  This 

2 

suggests  an  iteration  process,  the  next  step  being  to  take  ( 2/tt )  Yq  =  Y^  as  a  new 
interval  function  containing  -y",  which  leads  to  an  improved  YQ  and  a  corresponding 
value  for  the  interval  Rayleigh  quotient.  Indeed,  if  A^  C  AQ ,  then  the  exis¬ 
tence  of  an  eigenvalue  AG  of  A  is  guaranteed  by  the  Schauder  fixed  point  theo¬ 
rem  { 5  ]  . 

6.  Variational  problems  in  several  dimensions.  The  extension  of  Theorem  3.1 

to  problems  in  several  independent  variables  follows  immediately  from  the  corres- 

v 

ponding  extension  of  the  interval  integral.  In  R  ,  let  x  *  ,£2 , • • • r £  ) r  and 

the  region  of  integration  be  denoted  by  Q.  Following  the  prescription  given  in 

(1],  partition  (2  by  elements  Q.  ,  (2  , . , . ,  (2  with  measures  (areas  or  volumes) 

12  m 

i  =  l,2,...,m,  and  let 

(6.1)  VY.  =  [  (Y(x)l  ,  (Y(x)H, 

i  xGft .  . 

l  l 
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where  Y  is  an  interval-valued  function  defined  on  ft. 
of  all  partitions  of  ft  into  m  subregions,  then 
m 


(6.2) 


Z 

m 


n  T  y. -dft. ,  m  =  1,2,3,..., 

P  i=l  1  1 

m 


If  P  denotes  the  set  of 
m 


form  a  nested  sequence  of  closed  intervals,  and  thus  the  interval  integral  of 
Y  over  ft, 

(6.3)  L Y(x)dfl  *  n  £ 

ft  „  m 

m=l 

exists  for  arbitrary  Y.  It  is  not  difficult  to  show  that  this  interval  integral 
is  an  inclusion  monotone  function  of  its  integrand,  using  the  same  arguments  as 
in  11). 

Now,  one  can  let  denote  the  vector  of  partial  differential  operators  of 
order  i  in  RV,  for  example,  D  *  Q/3£  ,  3/3£  , .  .  .  , 3/3£  )  and  consider  the  func- 

1  12  V 

tional 


(6.4)  f  [y]  =  Jfif(x,  y.D^.D^ . D^dfl, 

v 

which  is  the  analogue  of  (1.6)  in  R  .  If  F  is  an  interval  inclusion  of  the  inte¬ 
grand  of  (6.4),  and  X  =  f(y]  is  a  stationary  value  of  f,  then  interval  integration 
provides  the  following  result. 

Theorem  6.1.  If  £  is  a  stationary  point  of  f  and  interval  vector  functions 

Y  ,  Y  ,  . . . ,  Y  exist  such  that  D.?e  Y.,  i  *  0,1, . . . , n  on  ft ,  then 
0  1  n  ii 

(6.5)  \  =  f [y ]  e  JSJF(x,Y0(x),Y1<x),...,Yn<x))dn. 

As  an  application  of  this  theorem,  suppose  that  in  the  values  of  y  are 
prescribed  on  the  boundary  3ft  of  a  region  ft,  and  one  wishes  interval  bounds  for 

(6.6)  X  =  min  } +  32y/3C2)dft 
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over  some  class  A  of  admissible  functions.  A  construction  similar  to  the  one 
in  §3  can  be  used,  or  Y^  can  be  constructed  on  the  basis  of  an  approximate  solu¬ 
tion  of  the  Euler  equation  for  (6.6) ,  which  in  this  case  is  simply  the  Laplace 
equation 

(6.7)  Ay  «  0,  y  -  yQ  on  30. 

I 
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